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We obtain a set of general formulae for determining magnetizations, including the usual electro- 
magnetic magnetization as well as the gravitomagnetic energy magnetization. The magnetization 
corrections to the thermal transport coefficients are explicitly demonstrated. Our theory provides a 
systematic approach for properly evaluating the thermal transport coefficients of magnetic systems, 
eliminating the unphysical divergence from the direct application of the Kubo formula. For a non- 
interacting anomalous Hall system, the corrected thermal Hall conductivity obeys the Wiedemann- 
Franz law. 

PACS numbers: 72.10.Bg, 75.10.-b, 75.47.-m, 66.10.cd, 44.10.+i 



Thermal Hall effect, or Righi-Leduc effect, is the ther- 
mal analogue of the Hall effect It gives rise to a 
transverse heat flow when a temperature gradient is ap- 
plied. There are recent surging experimental interests in 
studying the thermal Hall effect in various systems, re- 
vealing such as the phonon Hall effect [5, 's'l , magnon Hall 
effect [4], and so on. There are also many theoretical ef- 
forts to study these phenomena However, most of 
these theoretical studies face a fundamental issue: direct 
application of the Kubo formula, when done correctly 
without questionable "tricks" [2^ , alwaysyields unphys- 
ical divergence at the zero temperature Such an 
issue is actually a major obstacle in the theoretical stud- 
ies of the thermal Hall effect. 

The underl ying reason of the issue had been previously 
identified [III ll2| : in a system breaking the time-reversal 
symmetry, either by applying an external magnetic field 
or due to the spontaneous magnetization, the tempera- 
ture gradient not only drives the transport (heat) cur- 
rent, but also drives the circulating (heat) current that 
is not observable in the transport experiment. Both con- 
tributions are present in the microscopic current density 
calculated by the standard linear response theory, and a 
proper subtraction of non-observable circulating compo- 
nent is necessary. For the electric transport, such sub- 
traction involves the electromagnetic orbital magnetiza- 
tion density, while the subtraction of the energy current 
will involve the garvitomagnetic energy magnetization 
density which characterizes the circulating energy 
flow. However, the previous theoretical discussions do 
not clarify what the transport current and the magne- 
tizations are, and how the magnetizations can be evalu- 
ated for a general extended system. The issue becomes 
more fundamental because the magnetizations are gauge- 
dependent quantities |14] . and it is not a-priori clear 
what the proper gauges of the magnetizations should be 
when calculating the transport coefficients . 

In this Letter, we attempt to build the theory of ther- 
mal transport of magnetic systems on a firmer basis. 
We obtain a set of general formulae for determining the 



magnetizations, including the usual electromagnetic or- 
bital magnetization as well as the gravitomagnetic energy 
magnetization [Eqs. ((7HT0|]. We further show that these 
magnetizations naturally emerge as corrections to the 
thermal transport coefficients, recovering the Onsager re- 
lations and Einstein relations [Eq. (PO)) ]. and eliminating 
the unphysical divergence. The result is a complete set 
of general formulae for calculating the transport thermal 
Hall conductivity, as well as the other thermal-electric 
responses such as Nernst effect and Ettingshausen ef- 
fect [lit . The formula also clarify what the gravitomag- 
netic energy magnetization is and how it can be calcu- 
lated, and its thermodynamics is determined. We test 
our theory by calculating the thermal Hall coefficient of 
a non-interacting anomalous Hall system, and observe 
the emergence of Wiedemann-Franz law, consistent to 
the recent experimental observation (l6| . 

Preliminaries: To make our discussion specific, we con- 
sider a general electronic system. We should note that 
the formulae we will develop are general, applicable to 
the other systems such as the phonon and spin systems. 

We assume that the total Hamiltonian of the unper- 
turbed system can be written a.s H ^ J drh{r), where 
r denotes the spatial coordinate, and h{r) is the lo- 
cal energy density operator. To study the electric and 
thermal responses, we introduce the external mechanic 
fields: the potential (j){r) and the gravitational field ^/'(r), 
where the gravitational field is introduced as the me- 
chanic counterpart of the temperature gradient, following 
Luttinger fl3|. In the presence of these fields, the local 
energy density operator of the system is modified to [l3| : 



V^(r) = [1 + V(t-)] h{r) + (j){r)n{r) 



(1) 



where h{r) is the local density operator, and the Hamil- 
tonian of the system is -ff^,^ = / drh^^^(r). 

The particle and energy current operators of the sys- 
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tern are defined by the conservation equations |12l |: 

^ ^ ^[n{r),H,,^] = -V ■ J'^/'ir) , (2) 

^ yi^-'Hrl H,^„] = -V ■ Jt^Hr) , (3) 

where jf^"^ and J^''^ are particle and energy current op- 
erators, respectively. 

We further require that the current operators in the 
presence of the external fields can be related to the zero- 
field current operators Jjv and Je by [l^ : 



where 



jti^{r) = [l + ^,{r)]JN{r)., 



(r) = [l + V(rr JE{r) + <l>{r)JM{r) 



We note that the current operator is only defined up 
to a curl by Eqs. ([IHS]). As we will show later [See 
Eq. ()22p ]. one may use this freedom to find appropriate 
forms of current operators that do satisfy these scaling 
relations [26|. 

When the system is in equilibrium and in the absence 
of the external fields, we have V- 

•^N{E)i''') — 

particle (energy) current for the equilibrium density ma- 
trix /5o = (l/^o) cxp 



V- J^'' = 0, where 
N{E)if)') is the expectation value of the 



-K/UbTo , where K = ^ drK{r) 

and K{r) = h{r) — fion{r). As a result, we can introduce 
the zero field particle magnetization density A4'jv(r) and 
the energy magnetization density Msir) so that: 



The equation can also be considered as the (incom- 
plete) definitions of the magnetizations. To make the so- 
defined magnetizations physically meaningful, one needs 
to further require the magnetizations being the proper- 
ties of material, i.e., they should be well-behaved func- 
tions of r, and vanish outside of the sample. We also 
introduce the zero-field heat magnetization: MQ{r) = 
Msir)- noMNir). 

Magnetizations: We rigorously prove that, with the ap- 
propriate current operators that follow the scaling laws 
Eqs. (IIHS]), the total magnetizations can be calculated 
from the following set of equations: 



Mat -To 



2Mq-To- 



dfio 
dTo 
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^ Vq X (X_q; Jw,q 



— Vq X {fl_g;jQ^q 



— Vq X {K_q-jQ^q 



q^O 



(7) 
(8) 
(9) 
(10) 



(l//?o) /(f" dATr paa{-ihX)b 



the 



Kubo canonical correlation function (/3o = l/fcsTo) \li 



^N(Q) = JdrMN(Q){r), Jgir) = Jsir) - fioJN{r), 
and fiq, Kq, JN,q, JQ,q ^'^'C the Fourier transform of n(r), 
K{r), JN{r), Jgir) [dq = J dra(r)e^"' ''), respectively. 

Equations (|7l-fT0| are the central results of this Letter. 
The total magnetizations can be obtained by integrating 
over either the chemical potential fiQ [Eqs. ([71 [!])[ or the 
temperature To [Eqs. (|51 ll0p [. The corresponding bound- 
ary conditions are that at fiQ -^■ — oo, M^i^q-^ — and 
at To — > 0, {2Mq) coincides with right hand side 

(RHS) of Eq. ©[(Iinill, respectively. For electronic sys- 
^4~j tem, the two approaches are equivalent. On the other 
hand, for systems without the chemical potential, such 
(5) as the phonon and magnon systems. Equation ([TO)) is the 
only option for calculating the heat (energy) magnetiza- 
tion. 

In Ref. [l^, a similar formula for the electromagnetic 
orbital magnetization M = ^cMn was derived from its 



thermodynamic definition M — —{dil,/dB)Tg,fj.o, where 
il is the grand thermodynamic potential, and B is the 
magnetic field. It is easy to identify that RHS of Eq. ([5]) is 
nothing but -{dK/dB)^„^To , where K = n + TqS and S 
is the entropy of the system. Similarly, Eq. ([7]) is just the 
Maxwell relation between dM/dn and dN/dB, where N 
is the total particle number of the system. 

One can develop a similar thermodynamic interpreta- 
tion for the heat magnetization as well. For this pur- 
pose, it is necessary to introduce a fictitious "magnetic 
field" Bs which couples to Ms = Mq /Tq so that Ms = 
~{dQ./dBs)^_io,To- Bs can be related to the physical grav- 
(6) itomagnetic field Bg [ISl by = -{To/c^)Bg. In anal- 



ogy to the particle magnetization, RHSs of Eqs. ([OHI 
are -Tf){dN / dBs)To,t^o and ~To{dK/dBs)To,tJ.o: respec- 
tively, and these equations are nothing but the thermo- 
dynamic relations. It is important to note that the par- 
ticular way to introduce the thermodynamic quantities 
(e.g., Ms instead of Mq) is necessary for accounting for 
the extra factor of 2 in front of Mq in Eq. PU)) . 

We sketch the proof of Eqs. ([T HIO)) in the following [21|- 
We introduce the static response functions: 



X,,(r,r') = /3o(A%(r');AJ,(r))^, z,j = 1,2, (11) 

where ni(r) = n(r), h2{r) = K{r), Ji(r) = Jjv(^), 
J2(r) = Jqir), and Aa = a - (a)^. Applying Eqs. 
we obtain V-Xij(r,r') = (1/ifi) {[njif'), niir)])^, which 
implies: 

V • xUr) + iq ■ [xUr) - V x M„ (r)] = , (12) 

where xfjir) = J dr'xtj{r ,r')e-'i- \ Mn{r) = 0, 
Mi2(r) = Mjv(r), M2i(r) = MAr(r), M22(r) = 
2MQ{r). In deriving Eq. (HH), we have utilized Eqs. ^ 
13]) and Eqs. (|3H51)i which imply the operator form of the 
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commutators [fij(r'), fii(r)]. Equation ^ is then used 
to determine the equihbrium expectation values of the 
resulting commutators. 

Therefore, Xiji^) rnust have the decomposition: 

xUr) = iq X M,,(r) + e-'i'^V x K^,(r) . (13) 



due to the inhomogeneous local chemical potential and 
temperature field. We assume that the deviation from 
the homogeneity is small so that /i(r) w /io + Sfi{r), 
l/r(r) w (l/To) + d[l/T{r)]. By applying the static 
response theory we obtain, to the linear order of 
xi(r) = 5fi{r) and X2ir) = -Tod[l/T{r)]: 



Because both Mij{r) and Xiji^) ^I'e properties of ma- 
terial. nfj{r) must also be well-behaved, i.e., it should 
be bounded and vanish outside of the sample. Moreover, 
the value of ^^^(r) at the long wave limit (q — 0) can 
be related to the macroscopic thermodynamic quan tities 
[aM^(Q)(r)/aMo]j,„ and [aM^(Q)(r)/aTo]^^ M- Ap- 
plying (i/2)VqX to both sides of Eq. ([T3| . taking limit 
q — > 0, substituting Mij and i^l^'^ and integrating over 
r, we obtain Eqs. (ff HTO)) . 
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E 

3=1 



dr'x^,(r-,r')x,(r'), (17) 



where Xij(T', r') is the static response function defined in 
Eq. (HH), and jp{r) = ^jf''''{r)'j , which can be deter- 
mined by Eq. dH) and Eqs. (HHS]). Substituting Eq. ^ 
into Eq. ([T7]), and after some algebra, we obtain, to the 
linear order of (/), ^, 5^ and 5{\/T) [2l|, 



Thermal transport coefficients: We can show that the 
magnetizations determined by Eqs. ([7HTU1) will emerge 
naturally as corrections to the thermal transport coeffi- 
cients. To see this, we calculate the full response of the 
currents to small deviation from the global equilibrium. 
In this case, the system can be approximately described 
by the density matrix: 



P ~ P\cq + Pi 



(14) 



where pieq is the local equilibrium density matrix char- 
acterized by the local chemical potential /i(r) and local 
temperature T{r): 



1 

Picq = cxp 



- / dr 



h{r) — ^{r)h{r) 



kBT{r) 



(15) 



pi is the linear response correction to the local equilib- 
rium density matrix, determined by the Liouville equa- 
tion ihdp/dt + [p,H4,.^] = We define a(r) = 
[1 + iir)Mr) + fi{rjl P{r) = l/fcB[l + ^(r)]r(r). It 
is easy to see that when a(r) and /3(r) are spatially 
uniform, pieq becomes the exact global equilibrium den- 
sity matrix corresponding to the Hamiltonian i?^,^, and 
pi = 0. Therefore, the conditions of the global equilib- 
rium are Va{r) = and V/3(r) = 



We define jf'"^ = J^'^ and X 



N 



a{r)jfj^' . The forces conjugate to these currents are 
Xi = —(3{r)'Va{r) and X2 = V/3(r), respectively, so 
that the entropy generation is ds/dt + 'V ■ {(SJq^) = 

J2i jf'^ ■ [l]- The expectation values of the currents 
have two parts of contributions: 



floq 



rKubo 



(16) 



where J^^^° = Tipijf'^ is just the usual linear re- 
sponse contribution with the response coefficients deter- 
minable by the Kubo formula |2^. Besides this, there 



is an extra contribution J; 



Icq 



Trpioq jf ''^ , which is 



1 



4°^(r) «V X M'^-^'ir) - a{r)V x Mp^{r 



(18) 



^MN{r) xXi~ |MQ(r) x X2 , (19) 



where Mf/^{r) 



[l + ij{r)]MN{r) + SMN{r), 
1 + V(r)]^ [MEir)+(t)ir)MNir)]+SME{r), 
and SM]y^E){'i') is the correction to the particle (energy) 
magnetization due to the spatial gradients of the chemi- 
cal potential and temperature, determinable by Hifj{r). 

Applying Eqs. ([TCI [T51 [TO)) . we can obtain the total cur- 
rents responding to the non-equilibrium forces. However, 
due to the presence of J,-'^'^, such responses break the fun- 
damental non-equilibrium thermodynamic relations [l[: 
(1) Onsager reciprocal relations; (2) Einstein relations, i. 
e., the currents should only be proportional to Va and 
V/3, and vanish when the system is in the global equi- 
librium. The problem can be remedied by defining the 
transport currents as jf/^E) 



T<t>,-4> 



-VxMf^fj^y and 



the corresponding transport responses then become: 







.J2. 





^(11) 

^(21) _ Miv X 



^(12) _ JVf„ 
7^(22) _ 2A4-Q 



Xi 

X2 



(20) 

where Jf = (l/V) / drjf -''''^'{r), and V is the total 



- m 



volume of the system. L ^'■^^ is a tensor of rank two 



with the component Vai^' — dte (^Jj.-y; ji.a{t)'^ 

(a, 7 = x,y,z), which is the usual response coefficient 
determined by the Kubo formula [22j. It is easy to verify 
that both the Onsager relations and the Einstein rela- 
tions are recovered. The magnetizations determined in 
Eqs. ((T HTO)) naturally emerge as the corrections to the 
thermal transport coefficients. 

Application: We can apply these general results to study 
the thermal Hall coefficient of a non-interacting anoma- 
lous Hall system [isl . [23| , and show how the unphysical 
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divergence is eliminated and the Wiedemann-Franz law 
emerges. The energy density of such a system, in the 
presence of the external fields 4>{r) and ip(r), can in gen- 
eral be written as: 

+ip^{r)[V{r) + (t>{r)]cp{r)} , (21) 

where (p{r) {(f^{r)) is the electron annihilation (cre- 
ation) field operator with the two spin components, 
V = (l/m)[— iftV + Aso{r)] is the velocity operator with 
Aso{r) being the non-abelian gauge potential charac- 
terizing the spin-orbit coupling, and V{r) is the pe- 
riodic potential. The field operator ip{r) satisfies the 
Schrodinger equation: ihdip/dt — H^^iptf with Ti.^^^ = 
{m/2)v ■ [1 + ipir)] v + [l + ip{r)] [V{r) + 0(r)]. An ap- 
propriate energy current operator that does satisfy both 
Eq. (UJ and the scaling law Eq. ^ is, 



(r) 



1 + 



8^^ 



(1 + V)^(t)(^)^ X (i)^) . (22) 



The presence of the last term is essential for satisfying 
the scaling law Eq. 

With the appropriate energy current operator at hand, 
we calculate the thermal Hall coefficient. The usual Kubo 
formula yields. 



,Kubo 



where 11; 



nk 



Im 



r(22) 



dUn 



2TohV 



(23) 



^nk 



(Tik + ^nk — 2^oj 

is the periodic part of Bloch wave function for band n and 
quasi- momentum k, fnk = fi^nk) is the Fermi distribu- 
tion function, iik = {l/2m)[-ihV+Aso{r) + hkY+V{r), 
and tnk is the electron dispersion Q. It is easy to see that 
the coefficient diverges at the zero temperature. 



We 

(/3o/2i) 



calculate 



M, 



Q 

. and obtain: 



nk 

+ 2f^:;,(e„fc-/io)'/;fc|, (24) 

where = —2lm{dunk/dkx\dunk/dky). Mq is ob- 
tained by integrating Eq. (llOp . After some algebra, we 
(2Af§/Toy): 



obtain K% = A^Kubo 



"xy 



de(e-/xo)V,,(e)/'(e), (25) 



where (Txy{^) = Se„fc<e ^nfc ^^'^o temper- 

ature anomalous Hall coefficient for a system with the 
chemical potential e [isL [23| . It recovers the Wiedemann- 
Franz law at the low temperature fesTg ^ /^o 11, IH], 
and the unphysical divergence is eliminated. 

In summary, we have developed a systematic approach 
for calculating the particle and heat (energy) magnetiza- 
tions. We also explicitly show that these magnetizations 
naturally emerges as the corrections to the thermal trans- 
port coefficients, recovering the Onsager and Einstein re- 
lations, and eliminating the unphysical divergences. Our 
approach make no assumption on the nature of the sys- 
tem, so is equally applicable to fermionic (e.g., electron) 
or bosonic (e.g., phonon, magnon) systems, either non- 
interacting or interacting. The approach does not involve 
the ill-defined spatially extended operators, so is usable 
in practical calculations. 
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I. KUBO'S CANONICAL CORRELATION FUNCTION 



Kubo's canonical correlation function is defined as 



B-A) =^Tr 
Po 







(SI) 



where is the Hamiltonian of the system, and pcq ~ {^/Z) exp (^—(3qM'^ . Some of its properties used in the mair 
text are 



and: 



/3o ( AA; AS )^ - - /?o ^ AA; AB^^ - - 



AA, AS 



where: 



XBA = Poi^B-^A 



(S2) 



(S3) 



where Ai = i - ^i^ and A = {l/iK)[A, 

For a system perturbed by a static external force: 

^' - Mx. (S4) 

The change of the expectation value of an operator B can be calculated to the linear order [l| : 



(S5) 



(S6) 
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II. DETAILS OF DERIVATION FOR EQS. (7-10), MAGNETIZATION FORMULA 



The derivation of Eqs. (7-10) is detailed in the following: 

(1) We introduce Xij(r,r') = Pq (jS.nj [r') ; AJi{r)^ with i,i = 1,2, where fii(r) = n(r), h2{r) = K{r), 
Ji(r) = Jw(r), J2{r) = Jqir), Aa = d - (a)g and (■ • • )g = Tr [pa • • ■], so we have: 



V ■ X^J {r, r') =/3o (An, (r') ; V • A J, (r) ) , 



= - f3o {fij (r') ; -hi (r) 

From Eq. jSTI to Eq. jSSl, we have used V • Ji (r) = -iii and V ■ = 0. Using Eq. jSl]), we obtain: 

V • X»j (r, r') = ^ ([fij (r') , (r)])o . 
We define xlj = / dr'xij {f, f') e^"' ^''^'^ ), and have: 

V • X^, {r) + iq ■ Xl (^) - ^ / d^' d". (^') ' (^)])o , 



(S7) 
(S8) 

(S9) 
(SIO) 



where «, j — 1, 2. 

(2) We can obtain the right hand side of Eq. (ISlOp from the definitions of currents and their scaling laws. Basically, 
we have: 



^[n{r),H*'^]=-V-jtjHr)^ 



and: 



jr(r) = [l + VXr)]JA.(r), 

J^'^(r) = [1 + 4'{r)f \jE{r) + (r) J^(r)' 



If we set (r) = and 1 + -0 (r) = e"^ '', Eq. (jSlip becomes: 



e"J '^ /i [r') , fi (r)J = V • (^e'-J"- Jjv (r) j , 

= e'-'-'-iq • Jjv (r) + e'" '^ V • Jjv (r) 



and we obtain: 



1 



dr 



' (J [e''^-('''-'')/i (r') , fi (r)] = iq • V x (r) , 



(Sll) 
(S12) 

(S13) 
(S14) 

(S15) 
(S16) 

(S17) 



where we have used V • J'^ (r) — and J'^ (r) = V x Mn {r). This is exactly the right hand side of Eq. (jSlOP for 
z = l,j=2. ^ 

Using the similar approach, we can prove that: 



i I dr' {[n, (r') , (r)])^ e-^'^i^-'^') = iq ■ V x M,, (r) , 

where Mn{r) = 0, Mi2(r) ^ MN{r), M2i(r) = MAr(r), and M22(r) 2MQ(r). 
(3) Therefore, Xiji^) satisfies the equation: 



and it has the general solution: 



V ■ Xlj (r) + iq ■ [X^, (r) - V X M^- (r)] = , 



Xl, (r) = -ig X M,, (r) + e-^i-^V x k^^ (r) , 



(S18) 

(S19) 
(S20) 
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where K^Ar) is an arbitrary function. This equation can be considered as a decomposition of xfji'^)- It is important 
to note that the decomposition is not necessary to be unique, because the magnetization can only be defined up to a 
gradient. However, the arbitrariness does not affect our result on the total magnetizations, as long as both Afy (r) 
and Kf- [r] are well behaved functions: i.e., they are bounded for all r. The constraint is necessary because, first, 
magnetizations are properties of materials; second, only when these functions are well behaved, can their contributions 
presented in Eq. (18-19) be well defined. 

(4) We can relate (r) to the macroscopic thermodynamic quantities. To see this, we use Eq. (17) and see how 
the equilibrium currents are perturbed by the spatially uniform changes of the chemical potential and the temperature. 
We have: 



5Jr{r)^ I dr'[x.i(r,r')<5A*o-X.2(r,rW(l/To)] , 



V X 



K^,-°ir)6t,o-^^li''{r)To6{l/To) 



1=0 ( 



(S21) 

(S22) 
(S23) 



Note that SJ^ir) = Tr j|''^(r)5picq - a{r) jf^'"'' {r)Spicq « SJ^^^ir) - no5J^^{r). On the other hand, SJ^"^ is, by 
definition: 



Comparing the two sides, we obtain: 



Sjp = V X (SMn) , 

^ = V X (SMe) - /ioV X {6Mn) 



OMn (r) 



OMn (r) 



To 



dTo 
OMq (r) 



(.r) =To 



dMQ{r 



Mn (r) , 



dTn 



(S24) 
(S25) 



(S26) 
(S27) 
(S28) 
(S29) 



(5) Equation (|S20p can be rewritten as: 

Xl, (r) = iq X [M,, (r) - e-'''"-^ (r)] + V x [g-''^-'-^ (r) 
Applying VqX to the both sides of Eq. (|S30|) and setting q ^> 0, we obtain: 



(S30) 



2 V, xxf,(r)|^^g 



-M„- (r) + nl-" (r) - V x C7„- (r) 



(S31) 



where Uij (r) = ^Vq x (e "' '"/^fj (^•')lq^o' ^^^^-"^ substituting different components of Mij (r) and /^^^ (r) and 
integrating over r we come to the formulae for the total magnetizations: 



Mn-To 



2Mq - To 





_ 


dfio 


2i 


OMn 


_ /3o 


dTo 


2i 


OMq 




d^o 


2i 


dMg 




dTo 


2i 



-qi '^N.q 



Vq X (n_q; Jg^q 



Vq X (K_q;jQ^g 



q^O ' 
q^O 



(S32) 
(S33) 
(S34) 
(S35) 



In the derivation, we assume that / dr V x Uij (r) = 0. This is guaranteed because nf - is a well behaved function. 
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III. DETAILS OF DERIVATION FOR EQS. (18-19), LOCAL EQUILIBRIUM CURRENTS 

Inserting Eq. (13) into Eq. (17), we obtain: 



J^(r) « Jt\r) + ^ i^M,, (r) x Vx, (r) + J ^ V x k 

Because Xjq = / dr'xj (r') e^"' '' , we have: 

Jl'='^(r) « jpir) + ^ (^M,, (r) x Vx, (r) + V x | dr'^,, (r, r') x, (r' 



(r)x,J . (S36) 



=V X 



(S37) 



where {r,r') = J dq/ (2^)^ k^^ (r) e-'-J '"'. 

We can obtain J°'^ (r) through the scahng law. Without tp and 4>, we have: 

J7iE)i^)=^ ^M^iE){r). (S38) 

When -0 (r) and (p (r) are present, according to the scahng law in Eq. (|S13[) and (jS14[) we have: 

Ji'^(r) = [l + 7/>(r)] V X Mw(r), (S39) 
Ja" (r) = [1 + V(r-)]^ [V x Msir) - /i(r)V x MA,(r)] . (S40) 

For i = 1, inserting Eq. ([SSQ]) into Eq. ((S37| . 

2 

J}'^'i(r) « [1 + V(r)] V X MAr(r) - Mjv (r) x To + V x /" dr' ^ /^i, (r, r') x, (r') , (S41) 

2 

= Vx([l + VXr)]Mjv(r))-^MAr(r)xX2 + Vx /" dr' ^ /^i, (r, r') a;, (r') , (S42) 

so we can write: 

J'^^ir) « V X Mp^ir) - ^M^(r) x X2 , (S43) 

where M^^'^(r) = [1 + ^{r)] Mn (r. To, ^o) + 5Mjv(r) and ,5Mjv (r) = E'=i / ^^'/^y r') x, {r'). 
Similarly, for i = 2, inserting Eq. (jS40p into Eq. (|S37p . 

4'"^(r) « [1 + V(r)]^ [V X ME{r) - ^iir)V x MAr(r)] + Mat (r) x - 2Mq x TqV^ (S44) 
+ V X f dr'Y^ K2j (r, r') (r') , (S45) 
(1 + ^ ir)f (Mis(r) + (r) Mjv (r))] - a(r) V x ((1 + ^ (r)) Mjv(r)) (S46) 

2 

^M^ir) x Xi - |MQ(r) x X2 + V x ^ J dr'/«2j (r,r') (r') . (S47) 



Further, by substituting M'^''^ (r) into Eq. (|S46[1 we can write J2°'^ (r) as: 

4^'i(r) «V X Ml'^^lr) - a(r)V x M^^'^(r) - iM^(r) x Xi - ^Mgir) x X2 , (S48) 

where M^'^ {r) = (1 + (r))' (M^ (r, Tq, ^o) + <^ (r) Mjv (r, Tq, a^o)) + (r), -JM^ (r) 

/ c^''''^2j {'^^ '^') Xj {r'), and K'a^- = K2j + A*o«^ii- 
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IV. DETAILS OF DERIVATION FOR EQ. (22), DEFINITION OF ENERGY CURRENT 

The energy density can be written as: 

/i*^'^(r) = [1 + ?A(r)] |y [vifiir)]^ ■ [vip{r)] + ip^r) [Vir) + 0(r)] (^(r)} . 



(S49) 



The Schrodinger equation for the system is ih^ = , where Ti'^''^ = fv-[l + ip{r)] * + [1 + ipir)] [V{r) + 0(r)] . 

Therefore, we have: 



+<^^(r) [T/(r) + cp{r)] \n^^^^^{r) \ - \n*''^'^{r)\ ^ [V{r) + 4>{r)] ^(r) 
- V • I ^ [1 + i'{r)] nt)^(r)]^ [?^;'^^'''(^(r)] + [?i'^''^(^(r)] ^ [t)(^(r); 



so we can identify J"^'^ (r) as: 



Because: 



where = we obtain the fohowing scahng equation: 

Jt"^ (r) = [1 + 4,{r)f [Je (r) + (r) Jn (r)] + V (1 + V^(r))' x A (r) , 

where A(r) = ^ ('"'p)^ x {v^)- 

In order to satisfy the scahng law Eq. (5), we redefine the energy current operator as: 

J^'^^M ^ Jl^'^(r)- Vx ((l + ^(r)f A(r)) , 
JE{r) Je (r) - V x A(r) . 



(S50) 
(S51) 
(S52) 

(S53) 
(S54) 
(S55) 



(S56) 
(S57) 



This is exactly the energy current definition Eq. (22). It is straightforward to show that modified energy current 
operator satisfies the scaling law Eq. (5). 

The particle current operator is defined as usual. It automatically satisfies the corresponding scaling law Eq. (4). 



V. DETAILS OF DERIVATION FOR EQ. (23), KUBO CONTRIBUTION 

The thermal current operator Jq^ [r] is: 



Jqx (r) = 



(r) (r) + {v^(fi 



- ^ E ^7 (i^'x^ {r))^ ^iV (r) ~ (i-f'P (r))^ v.^ {r)) , (S58) 



where "f = x,y, z and we have set (j){r) = and ip (r) = 0. According to our definition, we have [2] k 
and: 



r (22) 

Kubo — xy 



fnk fn'k' 



{lpnk\ jQy IV'n'fe') O'n'k'lJQx li^nk) , 



(S59) 
(S60) 



6 



where ipnk is the Bloch wave function for band n and quasi-momentum fc, /„fc = /(e„fc) is the Fermi distribution 
function, and e„fe is the electron dispersion. According to our definition Eq. (jS58l) for Jqx, we have: 



(V'n'fe'l JQx \^nk) 



Ki^nk ) ri ^ 

TT: [{y^Vxi^n'k' I^V'nfc) + {Vx4>n'k' |V^W^?/'; 



2 8i 

7 

(V^W^Vn'fc' \Vx'^nk) - {Vflpn'k' |V^?)j;Vrifc)] ■ 

.Klpn'k' \Vxi^nk) = i^n'k - Mo) (Wn'fe| Wfex l^nfe) (^fcfc' , 



Note: 



with Vkx — dUk/d (hkx), and: 

(V^VxIpn'k' \Vy1pnk) + {Vx->Pn'k' IV-yU-yVnfc) = " (V-'n'fe' I W^^V^U-y {ipuk) + (V'n'fc' I ^'x V^W^ jV'nfc) , 

-0, 

and similarly, (V^u^Vn'fe' I'SxV'nfc) + {vjfpn'k' lyjVxipnk) = 0, so we come to: 



^(22) 



E 
E 



(/nfc — /n'fc) (Crifc + ^n'k — 
4i (e„fc — Cn'k)^ 

fnk (Crifc + En'fe " 



(Unfcl Wfca |Un'fc) ii'n'klvkx \lpnk) , 



Im [(Wnfel Vky \Un'k) {Un'k\vkx \Unk} 



Using the identity: 



we have: 



(^^n'fcl Vkx \Unk) — {^nk — ^n'k) ( Un'k 



du 



nk 



nk ^ 



ik + <^nk — 2/i( 



dkx I ' 



dk,, 



Ink ■ 



The formula can be rewritten as the alternative form. We introduce the new notations: 



m2 



mi 



(^)-^i-E\^ 

(e)^-Im^ 



(e) 



Therefore, we can express K^y^° as: 



,Kubo 



2. / dunk 

fii-E\^ 



Hk-e 
dunk 



nk 



dky 
dUnk 



dk. 



^ (e - fnfe) 



dky 



S{e- e„fc) 



2roV^ 



m2 (e) + 4 (e - ^o) mi (e) - 2 (e - ^o) (e) / (e) 



It is easy to see K^y^° is divergent in the low temperature limit. 



(S61) 
(S62) 

(S63) 



(S64) 
(S65) 



(S66) 
(S67) 

(S68) 
(S69) 

(S70) 
(S71) 
(S72) 

(S73) 



VI. DETAILS OF DERIVATION FOR EQ. (24), ENERGY MAGNETIZATION 



To calculate Mq^z, we use 2Mq — 7o-gjr^ 
Ma 



We can show: 



d 



(S74) 
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So we have: 



Mr 



TT;— > '\Wnk\ -Z \WrL'k') \Wn'k'\jQx,q\WrLk) 

^ nk.n'K' 



We have a careful calculation of (V'n'fc'l jQx,q IV'nfc); 



(V'n'fe' I Jqi,<j IV'nfe) =" 



8i 



In Eq. (|S76p . we can show: 



K-ipn'k' 



(V^W-yl/'n'fc'l e \Vx-ipnk) - {v-tlpn'k'le {S/jV^lprik)] 



= {Un'k-q \ Kk-qVkx \Unk) Sk',k-q ■ 



In Eq. (|S77p . similarly: 



and: 



= ~ {Un'k-q\ Vk-qx (V^ + Ifc^ - i^^) Vk-y |u„fe) 6k' ,k~q ■ 



{Vx1pn'k'\e "^"^ IV-yW^V'nfe) = ( (V-y + i/c-y) -Cfc^ |u„fc) 5k',k-q , 

(V^-O^Vn'fc'l e""'''' |Wa;'/'nfe) = - (Wn'fe-ql i>k-q-y (V-, + ifc^ - IQ-f) Vkx \Unk) 5k' M~q , 
{Vjijn'k'l e""'''' IVjVx'ipnk) = (Un'fc-q| Wfe-q7 (V^ + Ifc^,) Vkx \Unk) Sk' .k-q ■ 



Therefore, {Tp„'k'\JQxe "'"^ \ipnk) is, 

{ipn'k'l JQxe'"'"^ \^pnk) = {Un'k-q\ 

Mq^z can be simplified as, 



Kk-qVkx + Vk-qxKk %7 ivk-qxVkj - Vk-q^Vkx) 



(S75) 



(S76) 
(S77) 
(S78) 

(S79) 
(S80) 



(S81) 
(S82) 



(S83) 
(S84) 
(S85) 



\Unk)5k'M~q- (S86) 



Mr 



d \ ^ fnk — fn'k~q , , Kk + Kk-q , , , , Kk — q^kx H" ^k—qx^k 



'idqy ^ 



[Unkl 



, . ^nk ^n'k~q 
nk.n' ^ 

T,-y fill i'>^k — qx'^k-y ^k — q'j^kx , 



2 

\Unk) ■ 
o 

First, we calculate AfQ*°' for n ^ n' . When q ^> 0, we have: 



l^n'k—q) {'^n'k—ql 



M'qf^-l i^nk + en'k-2fioflm 



7i^n'k 



Unk 



dUn'k 



dky 



{Un'k\ Vkx \Unk) 



fnk fn'k 
^nk ^n'k 



Using the identity Eq. (jS68l) . we finally come to: 



M inter 



nk 



du 



nk 



l7ik 



(S87) 
(S88) 

(S89) 
(S90) 



Next, we calculate Mq^^^ for n = n' . When q ^ 0, we have: 
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M '^T - 4 



nfc 



nk 



dUnk 



dk,, 



{Unk\vkx \Unk) 



dk 



KkVkx + VkxKk \Unk) 



f nk 
f 

J nk 



rife I VkyVkx l^rifc)] 



Using the identity Eq. (|S68p and after some simple algebra, we obtain: 



Ah ^ 

Therefore, we have: 



;\'>intra 



nk 



dUnk 



dkr 



{^^nk — iikj — 4 (e„fe — ^o) i^nk — 'Hk 



dUnk 



dk,, 



{^nk - Aio) fnl. 



2h 



Ah- 



^Im 


7 du„k 
\ dk^ 


(iik + ^nk 


nk 






^Im 

nk 


' 1 dUnk 

\ dkx 


^Cnfc — Hk 



dUnk 



dky 



fnk 



-Tik^ - 4 (e„fc - /^o) (cnfc - Hk 



dUnk 



dk,, 



(e,ifc - Ho) f'„ 



nk ' 



(S91) 
(S92) 
(S93) 

(S94) 

(S95) 
(S96) 



We use 2Mq_2 - To{dMQ, JdTo) = Mq^z to obtain Mq,^. Using the notations of Eqs. (|S70p - (|S72|) . we obtain: 

1 



Mr 



de 



-m2 (e) / (e) + 2 (e - ^^o) ™i (e) / (e) - 21^. (e) / dxx/ (x) 



(S97) 
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